Imaging quasi-particle wavefunctions in quantum dots via tunneling spectroscopy 
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We show that in quantum dots the physical quantities probed by local tunneling spectroscopies, 
namely the quasi-particle wavefunctions of interacting electrons, can considerably deviate from 
their single-particle counterparts as an effect of Coulomb correlation. From the exact solution of 
the few-particle Hamiltonian for prototype dots, we find that such deviations are crucial to predict 
wavefunction images at low electron densities or high magnetic fields. 
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Current single-electron tunneling spectroscopies in 
semiconductor quantum dots (QDs) [3 13- may pro- 
vide spectacular images of the QD wavefunctions, both 
in real 000] an d reciprocal space 0] • The mea- 
sured intensities have been generally attributed to the 
probability densities of ground or excited single electron 
states occupying the dot. As pointed out by Wibbel- 
hoff et al. |9|, however, the role of other electrons filling 
the dot may actually be relevant. Indeed, QDs can be 
strongly interacting objects with a completely discrete 
energy spectrum, which in turn depends on the number of 
electrons, A 0,0] • Therefore, orbitals can be ill-defined, 
losing their meaning due to interaction. Also, it is un- 
clear how many electrons one should take into account to 
calculate the total density of states, as a particle tunnels 
into a QD filled with A electrons. In this Letter we thus 
address the following basic questions: What are the phys- 
ical quantities that are actually probed by scanning tun- 
neling microscopies (STM) 000] or magneto-tunneling 
spectroscopies 0] of QDs? How do they depend on 
interactions? Can they deviate from the common single- 
particle picture in physically relevant regimes? If only 
one many-body state is probed at time, then the signal 
is proportional to the probability density of the quasi- 
particle (QP) being injected into the interacting QD. We 
demonstrate that the QP density dramatically depends 
on the strength of correlation inside the dot, and predict 
the wavefunction mapping to be a useful experimental 
tool to image QPs, both in direct and reciprocal space. 

The imaging experiments, in their essence, measure 
quantities directly proportional to the probability for 
transfer of an electron through a barrier, from an emitter, 
where electrons fill in a Fermi sea, to a dot, with com- 
pletely discrete energy spectrum. In multi-terminal se- 
tups one can neglect the role of electrodes other than the 
emitter, to a first approximation. The measured quan- 
tity can be the current 00, the differential conductance 
0E10III3 or the QD capacitance 0,0, while the emit- 
ter can be the STM tip 0, El , or a n-doped GaAs con- 
tact 000 E El , and the barrier can be the vacuum 
01| as well as a AlGaAs spacer H EH . 

According to the seminal paper by Bardeen the 



transition probability (at zero temperature) is given by 
the expression (2-k/K) \M\ n(e/), where M. is the matrix 
element and n(e/) is the energy density of the final QD 
states. The common wisdom would predict the probabil- 
ity to be proportional to the total density of QD states at 
the resonant tunneling energy, e / , possibly space-resolved 
since M would depend on the resonant QD orbital 
To proceed, let us assume that electrons from the emit- 
ter access through the barrier a single QD at a sharp 
resonant energy, corresponding to a unique, well defined 
many-body QD state, and reconsider the transition ma- 
trix element Mk,N for transfer of an electron from emit- 
ter to QD. Mk.N is given by [recasting Eqs. (6-7) of 
Ref. ^3 m second quantized form] : 



M k , N ({k},N -l\M\{k*},N), 



M = 



2m* 



.5* <9#t 



dz 



<5(zbar - 2) dr. (1) 



Here \{k}, A — 1) and \{k*},N) are two many-particle 
states of the entire system of similar energies, with A — 1 
and A interacting electrons in the QD, respectively, and 
the remaining A to t — A + 1 and A to t — A electrons, re- 
spectively, in the emitter. The fixed coordinate along the 
tunneling direction z appearing in l|TJl. Zbar, can be every- 
where in the barrier, and dr is the infinitesimal volume 
element. The fermionic field operator \&(r), destroying 
an electron at position r = (g, z), can be expanded over 
the basis of emitter and QD single-particle states, 4>k and 
4> a , respectively [l4|: ^( r ) = J2i ( f ) i( r )cii where i = k,a 
and we take unitary volume normalization. We omit spin 
indexes and summations for the sake of simplicity. We as- 
sume that electrons in the emitter do not interact and are 
associated to the two sets of quantum numbers {k} and 
{k*}, respectively, which differ in the occurrence of the 
index k labeling the electron which leaves the emitter and 
tunnels to the QD as |{fc}, A — 1) evolves to {fc*}, A). 
Moreover, we assume for convenience that the xy and z 
motions of electrons are separable, and that electrons in 
the QD all occupy the same confined single-particle state 
along z, XQd(z), namely (f> a (r) = Q (e)XQD(z)- Under 
these conditions we may factorize the matrix element as 
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Mk.N oc T k M k M, with 



T fc = 



2m* 



XfcO) 



dXQp(z) 
dz 



dz 



(2) 

where Xfc(z) is the emitter state along z evanescent in 
the barrier, <j>h{r) = 4>k{o)Xk{z), and 

Mk, N =J2[<t>Ue)<t>de)d e ({k}, N - 1 1 c{c a \{k*},N). 

Eventually assuming that the many-body states can be 
factorized into an emitter and a QD part, we obtain 



M, 



k,N 



(3) 



where (Pqb(q) is the quasi-particle (QP) wavefunction of 
the interacting QD system [lij: 



<Pqd(q) = {N-1\*(q)\N). 



(4) 




Results (HI II) are the key for predicting wavefunction 
images both in real and reciprocal space. In STM, (f>k(g) 
is the localized tip wavefunction; if we ideally assume it 
point-like and located at g [l3|. i.e. 4>k{o) ~ 5{g — Qq), 
then the signal intensity is proportional to </9qd( 
which is the usual result of the one-electron theory [ 
provided the ill-defined QD orbital is replaced by the 
QP wavefunction unambiguously defined by Eq. 
In magneto-tunneling spectroscopy, the emitter in-plane 
wavefunction is a plane wave, <f>k (q) = e ifc, e, and the 
matrix element Q is the Fourier transform of <^qd, 
Mk,N = <^QD(fc)- Again, we generalize the standard one- 
electron result by substituting Lpqi)(k) for the QD orbital 
[then Eqs. |gE) coincide with (A1-A2) of Ref.Q. Note 
that Mk,N is the relevant quantity also for intensities 
in space-integrated spectroscopies probing the QD addi- 
tion energy spectrum ll ll| . Consistently, in the non- 
interacting case, </9qd(£?) reduces to the highest occupied 
one-electron orbital 4> a (g) 0|: in this limit an electron 
tunnels from the emitter to the orbital <fia{o) which res- 



The 



onates at the Fermi energy, with \N) — c* a \N — 1). 
latter regime probably corresponds to most of the exist- 
ing experimental evidence 0, H, El 01 El ■ However, it 
is interesting to analyze realistic scenarios that deviate 
from the one-electron picture. 

Therefore, we study <£>qd((?) in a paradigmatic in- 
teracting case, and consider a few electrons in a two- 
dimensional harmonic trap, which was proven to be an 
excellent model for different experimental setups Q . The 
QD effective-mass Hamiltonian is 



N 1 
2T = 5>o(<) + 5 £ 



i*7 



K\Q t - Qj 



(5) 



with 



Ho(i) 



2m* 



p-^A(g)] 2 + m*u; 2 g 2 /2. (6) 



Here k is the static relative dielectric constant of the 
host semiconductor, and A(g) is the vector potential 
(A = B x g/2) associated with a static and uniform 
magnetic field B along z, which reduces the cylindri- 
cal spatial symmetry group of the system from -Dooh, at 
B = 0, to Cooh, when B ^ 0, making it chiral. The 
QD wavefunction has an azimuthal quantum number m, 
<Pqd{q) = ( ^QD(£ l )e lm ^, which is fixed by the total angu- 
lar momenta M of \N) and \N— 1), m = M^—Mn-i 1 and 
can be expanded over the basis of Fock-Darwin (FD) or- 
bitals ipnm(o) eigenstates of the single-particle Hamil- 
tonian ©: <Pqd(q) = Y^=Q a nVnmi.Q), where n's are 
radial quantum numbers and a n coefficients to be deter- 
mined. We solve numerically the few-body problem of 
Eq. (5), for the ground state at different N's, by means 
of the configuration interaction (CI) method |l6| , where 
\N) is expanded in a series of Slater determinants built 
by filling in the FD orbitals with N electrons, and consis- 
tently with symmetry constraints . Then, we evaluate 
the matrix element (@J), and find the values of a n for a 
truncated FD basis set. 

There are two ways of artificially tuning the strength 
of Coulomb correlation in QDs: one is to dilute electron 
density, and the other is to turn on B. In both cases, 
at low enough densities or strong enough fields, electrons 
pass from a "liquid" phase, where low-energy motion is 
equally controlled by kinetic and Coulomb energy, to a 
"crystallized" phase, reminescent of the Wigner crystal 
in the bulk, where electrons are localized in space and 
arrange themselves in a geometrically ordered configura- 
tion such that electrostatic repulsion is minimized [3ll20| . 

We first consider reducing the density at B = 0. The 
typical QD lateral extension is given by the characteris- 
tic dot radius ^qd = (h/m*^) 1 / 2 , £qd being the mean 
square root of g on the FD lowest-energy level ip 00 . As we 
keep N fixed and increase ^qd, the Coulomb-to- kinetic 
energy ratio A = ^qd/ob ~ ft 2 «/(m*e 2 ) is the effeo 



tive Bohr radius of the dot] increases as well, driving 
the system into the "Wigner" regime |2l|. As a rough 
indication, consider that for A ~ 2 or lower the electronic 
ground state is liquid, while above A « 4 electrons form 
a "crystallized" phase fisj |. 

Figure ^ shows </?qd vs. g, as up to six electrons are 
successively injected into a "liquid" QD with a realis- 
tic density of A = 2 0. The QD filling sequence is 
well known |lfl ll ll| , in analogy with the Aufbau princi- 
ple of atomic physics: in the independent-electron pic- 
ture (A = 0, dashed lines), ^qd is the highest-energy 
occupied orbital which is filled by the electron added to 
the dot. Howevever, Coulomb correlation significantly 
spreads the wavefunction (solid lines) and moves the QP 
peak towards the QD edge. The spreading is caused by 
the increase of weights a n of high-energy FD orbitals, 
as interaction is turned on; nevertheless, the behavior of 
ipQD around g « is dictated by its angular dependence, 
<PQD(g) p' m ', while it decays like exp(—g 2 /2£q D ) as 
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FIG. 1: Quasi-particle wavefunction (solid line) vs. g for dif- 
ferent (N — 1) — > N transitions, with A = 2. The dashed line 
represents the non-interacting orbital (A = 0). The ground 
states are (M, S) = (0,1/2), (0,0), (1,1/2), (0,1), (1,1/2), 
(0,0), for N going from 1 to 6, respectively. S z = (S z = 1/2) 
if N is even (odd). The norm of the A = 2 wavefunction is 1, 
0.84, 0.84, 0.40, 0.37, 0.73, respectively. Insets: Total ground- 
state charge densities n(g) (arb. units) for N going from 1 (top 
left) to 6 (right bottom). The length unit is ^qd. 



g — > oo. The QP amplitude is strongly suppressed in the 
(N — 1) — > N tunneling processes involving the N = 4 
open-shell ground state, with respects to other additions 
(Fig. (IJ. This is a spin-blockade effect, since the total 
spin, S, is maximum at N = 4 (S = 1 according to 
Hund's rule [13 )j an d we assume that its ^-component is 
zero, S z — 0. Besides, the general trend is that the QP 
wavefunction norm, hence the integrated experimental 
signal, diminish as N and A increase (see also Fig. |5J. 

Note that the interpretation of tunneling spec- 
troscopy in terms of the total density, n(o) — 
(N\&(q)^(q)\N)/N, is inconsistent with our point of 
view, as it is seen by comparing QP wavefunctions of 
Fig. with total densities for the corresponding N- 
electron states (insets). While total densities and QP 
probabilites resemble each other up to the addition of 
the second electron, after the third electron tunnels into 
the dot they can be clearly discriminated in the labora- 
tory: QP probabilities have a strong angular dependence 
(hybridizing degenerate states with ±m) and a node at 
the QD center, while total densities are approximately 




radius p (Z QD ) 



FIG. 2: Top: Quasi-particle wavefunction vs. q for different 
values of A, as the 6-th electron tunnels into the QD. The 
wavefunction norm, for A going from 0.5 to 10, is 0.97, 0.73, 
0.48, 0.32, 0.22, 0.15, respectively. Bottom: Six-electron to- 
tal density n(g) vs. g. The ground states for N — 5,6 are 
(Af, S) = (1, 1/2), (0, 0), respectively, for all A. 



circular (exactly, for TV = 4, 6) and filled. 

As one reduces the density, the appearance of QP 
wavefunctions dramatically changes. In Fig. we study 
the injection of the 6th electron, as A goes from 0.5 up 
to 10. Note that A = 0.5, 2, 4 (equivalent to GaAs lat- 
eral confinement energies hcu = 47, 3.0, 0.74 meV, re- 
spectively), typically correspond to different experimen- 
tal QD devices, such as self-assembled 0, 122V vertical 
mesa etched [l(| , and 2DEG-depletion QDs [lj ■ A six- 
electron Wigner molecule forms for A > 4, with one elec- 
tron localized at the QD center, and the remaining five 
arranged on an outer ring, at the vertices of a regular 
pentagon 0, [U 0] • The crystallization is clearly seen 
in the bottom panel of Fig. where, as A increases, the 
total density develops one peak at q = 0, for the cen- 
tral electron, and another one, close to g = 3^qd, for 
the outer ring. Similarly, the five-electron molecule is a 
hollow pentagon |l8|, |23j . In the top panel of Fig. |21 we 
see that the QP wavefunction is strongly affected by elec- 
tron localization: while for A < 4 it somehow resembles 
the non-interacting FD orbital (n,m) — (0,-1), being 
spread uniformly across the dot, for A > 4 it develops a 
well formed peak close to the outer-ring position. The 
QP weight in the region inside the ring is strongly de- 
pleted, eventually appearing as a shoulder of the main 
peak. We conclude that, in the crystallized phase, the 
6th electron can only enter the external ring, with neg- 
ligible probability of being located into the center. For 
smaller N, we find that electrons just enter the outer 
ring, since the pertinent geometrical phases are hollow 
regular polygons [2o|. 

We now come to the effect of a strong magnetic field 
parallel to the tunneling direction z. As B increases, the 
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TABLE I: Absolute value of the modulation coefficient |ao| of 
the quasi-particle wavefunction v 3 Qd(&) = aoVom(£?), where 
m = (N — l)/u, for different (TV — 1) —> N tunneling processes 
and filling factors v. 



kinetic energy is quenched, Landau bands of almost de- 
generate FD levels being formed. M increases due to 
Coulomb repulsion, since the higher m, the outer the FD 
orbital [2(|. In correspondence of "magic" values of M, 
the ground state turns out to be particularly stable [2^ : 
this family of incompressible |27j states has been vari- 
oulsy regarded as reminescent of fractional quantum Hall 
effect (FQHE) states in two-dimensional electron layers 
, or as a collection of Wigner molecules [28| . 

In analogy with FQHE, it is convenient to introduce 
the filling factor v, defined as v = N(N — 1)/2M, and to 
consider only FD levels in the lowest Landau band and 
full spin-polarization, which turns out to be a reason- 
able approximation at high B In realistic situations, 
there are significant B-ranges where v is constant as N 
is changed UH^. At v = 1. the interacting states are 
maximum density droplets 3, 30J , namely incompressible 
disks of almost uniform density, \N) = Ilm=o %m\®)> an< ^ 
ipQD is simply the highest occupied FD state, ipoN-i, lo- 
cated at the edge of the dot, which is being filled by the 
tunneling electron, with a n = 5 n o- 



<PQn(g) = (pn=Q,m=N-l((>)- 



(7) 



Equation Q is a remarkable result: while the total elec- 
tron density is a uniform disk, the measured squared 
modulus of QP wavefunction will be an annulus of the 
same radius as the charge distribution. If v < 1, the 
wavefunction will be still proportional to the FD or- 
bital f n =o.m, with m = (N — l)/u and ao =/= 1, namely 
</?Qd(£>) = ao </?n=o,m(£>)- The only effect of strong corre- 
lation in these regimes is to modulate the amplitude of 
the non- interacting wavefunction via the coefficient ao- 
Table [I] shows the calculated values of ao for various tun- 
neling processes at particularly stable filling factors. Ex- 
cept for some cases, |arj| monotonously decreases as v 
diminishes or as N increases. E.g., at v = 1/5 |ao| is re- 
duced by two order of magnitudes with respect to v = 1 , 
when the 6th electron enters the dot. Table [I] shows that 
interaction enforces very effectively the orthogonality of 
incompressible states |3jJ, and therefore we expect that, 
as a high field component is applied parallel to z, tunnel- 
ing is strongly suppressed by the reduction of the matrix 



element M/ Cj at [Eq. ©]: a purely many-body mechanism, 
the single-particle matrix element Tk [Eq. being left 
unchanged by the field. 

In conclusion, we have shown that quasi-particle wave- 
functions of QDs are extremely sensitive to electron- 
electron correlation, and may differ from single-particle 
states in physically relevant cases. This result is of inter- 
est to predict the real- and reciprocal-space wavefunction 
images obtained by tunneling spectroscopies, as well as 
the intensities of addition spectra of QDs. Close com- 
parison with experiment is not yet possible in the case of 
Ref. [2j , where many dots are probed at once and the con- 
finement is too strong. Promising samples are also those 
of Refs. 0, ^3 > allowing for access to a single dot and 
full control of N. We hope that our results will stimulate 
further experiments. We believe that our findings will be 
important also for other strongly confined systems, like 
e.g. nanostructures at surfaces |32j . 
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